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1. INTRODUCTION 

In urban areas, the motorcycle can be seen as the primary mode of transportation among road users, 
especially in Asian countries such as Malaysia, China, and Taiwan. Taking Malaysia as an example, 
according to the Road Transport Department Malaysia, the number of registered motorcycles contributes to 
almost 45.98% of the total registered vehicles by the end of 2015 [1]. The number is expected to increase 
further in the coming years as the public sentiment, especially in the urban area, is hugely affected by the cost 
of travelling from and to the home daily [2]. Evidently, according to the latest report in June 2017, the 
number of motorcycles has increased to almost 6.93% from 2015 [3]. 

Globally, with the increasing number of motorcycles and other vehicles on the road, the 
transportation sector has become the main contributor to air pollution and worsening the environment by the 
increase in greenhouse gas emissions [4]. As the dominant vehicle in the transportation sector, the 
motorcycle is contributing to almost 10% of air pollutants in Taiwan [5]. As part of mitigation measures in 
environmental protection, the government aims to lessen the number of emissions coming from greenhouse 
gas (GHG) and air pollutants produced from vehicles [6]. It has been studied that lower noise pollution, 


Journal homepage: http://ijpeds.iaescore.com 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 Oo 31 


zero-emission of GHG and 75% reduction of CO2 emission can be achieved through electric motorcycles 
which consequently act as an effective environmental mitigation measure in the urban area [7], [8]. 

The electric vehicle is more appealing to the consumer in low-power usage as it is a more 
cost-effective and more energy-efficient vehicle [9]. The development of low-powered electric vehicles at the 
moment has been focusing on an electric motorcycle (EM) due to the pressure exerted from the energy crisis as 
well as the worsening environmental impact caused by the usage of motorcycles in an urban area [7]. Indeed, 
the researchers have shown that EM is highly efficient for energy consumption as compared to the internal 
combustion engine motorcycle (ICEM) with the energy efficiency of EMs is about 40% more than gasoline 
motorcycles [10]. Even though the electrifying the vehicle can significantly contribute to the reduction of air 
and noise pollution, it does not add to a lower number of the motorcycle on the road regardless of electric or the 
conventional motorcycles. The road fatality is found to be the highest in Malaysia as compared to other 
association of Southeast Asian Nations (ASEAN) countries with 50% fatalities in road accidents involving 
motorcyclists [10]. Numerous studies identify human factors as the cause of the accident [11]-[13]. 

Carelessness, distraction, inexperienced drivers as well as the failure of controlling speed are some 
of the human errors that contribute to a road accident [14], [15]. However, there are also other causes for a 
road accident in which relate to vehicle defects such as mechanical failures, tire, and brake defects as well as 
many different potential reasons [16], [17]. The road environment also can add the potential factors to cause 
the road accident by the presence of an adverse condition of weather, road, and traffic [18], [19]. Increasing 
road fatalities and injuries are mostly caused by the failure of motorcyclists in adjusting the speed of 
motorcycles during monsoon season on wet road conditions. It is indeed, for dry or wet bare roadway, an 
increase in accident rate is very significant as the tire friction reduces [20]. The incidences are caused by the 
high slip ratio on wet conditions as compared to the dry road. 

Slip is generally given as a percentage of the difference between the surface speed of the wheel 
compared to the speed between axis and road surface. According to [21], the factors of the relation between 
slip ratio and different road conditions have been disregarded in the few existing research on two-wheel 
electric motorcycle dynamics. The climatic condition in Malaysia in which experiences frequent monsoon 
and sunny days has been included as part of the case study in this paper. In the research conducted by [22], 
an accident prediction model is developed in the hope to reduce the motorcycle casualty on Malaysia’s 
primary road. The model is established to provide a proper and systematic road geometry as well as traffic 
census inventory for Malaysia’s road planning. On the other hand, a real test of [CEM has been conducted by 
the West Virginia University Institute of Technology to assess the performance of ICE motorcycles travelling 
on the dry road [23]. A similar assessment was conducted at the University of Lisbon on the Saldanha 
electric motorcycle to observe the speed of the vehicle in urban road conditions [24]. It is acknowledged that 
not all disasters are preventable but their impacts including all the underlying risks can be reduced. 

Hence, in this study, the models of the motorcycle are developed to investigate the maximum and 
suitable speed range of the motorcycle on the dry and wet roadway. Both ICEM and EM models of motorcycles 
are included to be under this study. The suitable speed range of the motorcycles is investigated as a means to 
ensure the stability of the motorcycle as the slip ratio varies according to road conditions. This paper consists of 
the modelling of motorcycles in section 2. In section 3, the simulation design is presented, and the discussion is 
presented in section 4. Finally, in section 5, the conclusion serves as a summary of the findings. 


2. MODELLING OF MOTORCYCLES 
2.1. Internal combustion engine motorcycle (ICEM) 

The development of an internal combustion engine motorcycle or ICEM requires important 
functional components such as generic engine, torque converter, inertia block, simple gear, nonlinear 
rotational damper, differential, vehicle body, and tires. The generic engine model acts as a reference to 
supply power to the ICEM. The torque of the engine can be generated from the throttle of the motorcycle. 
The equations for the generic engine are shown as (1)-(7). Next, the torque converter is connected to the 
generic engine. The torque converter has an impeller shaft as an input and turbine shaft as an output whereas 
the equations for the torque converter are shown as (8)-(12). 


P(O, Te) = Teg N) (1) 
nak (Ge) C2) © 
p(w) + pitw + pw? — paw? (3) 
Pik P= 73 = 2 (4) 
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pı + 2p2 — 3p; = 0 (5) 


tw = = (-p + (p2? + Ap,p;) 2) (6) 
was (7) 
R, =" (8) 
wg 
T 
aa (9) 
K = w//(T) (0 
Ry = RyRy (11) 
K = K(Ry) (12) 


The input and output of the torque converter are connected with the inertia block. The inertia block 
is used to ensure the motorcycle is accelerated in positive a direction. The equation for the inertia block is 
shown as (13). Then, the torque converter is connected to the simple gear. The gear ratio can be varied by 
using simple gear. The equation for simple gear is shown as (14). The simple gear equation shows the 
relationship between the number of teeth of the follower on gear and the number of teeth of the base on gear. 
Furthermore, the nonlinear rotational damper works as coupling, drives and damper. The equation for 
nonlinear rotational damper is shown as (15). Based on the equation, it shows a relationship between the 
damping torque and angular velocity. The output from the nonlinear rotational damper is connected with the 
differential. The differential will allow the driven shaft to move or to rotate with varying velocity. The 
equation for the differential is shown as (16). 


T=] (13) 
Gre = e/g = Ne/ (14) 
Tp = b\w + sign(w) - baw? + baw? + sign(w) - b,w* + bsw® (15) 
Wp = (1/2) gp (Wei + Ws2) (16) 


The proposed model for the vehicle body of the motorcycle is shown in Figure 1. The vehicle body 
is being modeled with longitudinal dynamics and motion. The equations for the vehicle body are shown as 
(17)-(22). By modeling the (17)-(22), the normal force acting on the tire and speed of the motorcycle can be 
identified. 


My = F — Fa — mg - sinf (17) 
Fe = (Fas + Fer) 18) 
Fa = CapAVy — W)? - sgn CV — W) (19) 
Fyr = eee on 
Fzr = eee Si 
Fzf + Fzr = mg ŽP 22) 


The study on vehicle dynamic behavior requires the presence of a tire model [25]. The magic 
formula tire model has been introduced by Pacejka, whereby this particular tire model is used for the analysis 
of tire behavior when driving on different road conditions since it is capable to deal with different slip 
conditions. The capability of the model in dealing with different slip conditions is majorly influenced by the 
retrieval of data on the tire force, tire moments, aligning torque, and longitudinal force. The equations for the 
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magic formula tire model are shown as (23). The coefficient used in magic formula is B, C, D, and E where 
B is stiffness factor, C is shape factor, D is peak value and E is curvature factor. The summary of symbols 


used in this modeling section is listed in Table 1. 


F, = f (k, F,) = F,D sin(C arctan(Bk — arctan(Bk))) (23) 
Figure 1. Proposed model vehicle body of motorcycle 
Table 1. Parameters of ICEM and EM systems 
Symbol Parameter of ICEM system Unit Symbol Parameter of EM system Unit 
V; longitudinal vehicle velocity km/hr Ron resistance on bridge Q 
T inertia torque Nm V, nominal voltage V 
T: engine torque Nm Vp back emf V 
Tp damping torque Nm k, back emf constant - 
P1P2P3 power demand coefficient - T torque Nm 
g) engine power demand coefficient - Te electromagnetic torque Nm 
Q engine speed rpm ki torque constant 
R, speed ratio - i motor current A 
wr turbine angular speed rad/s V battery voltage V 
wi impeller angular speed rad/s Vref reference voltage V 
Rr torque ratio - Vinin minimum reference voltage V 
K capacity factor - Vmak maximum reference voltage V 
J inertia kg.m? Vp back EMF Vv 
SFB follower-base gear ratio - k, back emf constant - 
Nr number of teeth of follower on gear -- a, B curve-fitting constant - 
Ng number of teeth of base on gear D duty cycle - 
bı, damping coefficient - Ho input voltage to H-bridge Vv 
bo, a b5 
wp angular velocity of longitudinal driveshaft rad/s Vpwm pulse width modulation V 
(PWM) voltage 
Wsi, WS2 angular velocity of the sun gears 1 and 2 rad/s Apwm amplitude of PWM signal - 
m vehicle mass kg y output current A 
Vw headwind speed km/hr Ti Tə torque on input and output Nm 
shaft 
h number of wheels on each axle - power on the input and output wW 
shaft 
Fy Fyr longitudinal forces on each wheel at the front and rear N gear ratio - 
ground contact points 
Fy, Fy normal load forces on each wheel at the front and rear N 
ground contact points 
b inclination angle degree 
a distance of front axles from the normal projection of m 
vehicle CG in the common axle 
b distance of rear axles from the normal projection of vehicle m 
CG in the common axle 
n number of wheels on each axle - 
A effective frontal vehicle cross-sectional area m? 
Ca aerodynamics drag coefficient - 
p mass density of air kg/m? 
B stiffness factor - 
Cc shape factor - 
D peak value - 
E curvature factor - 
F.=f(k,F,) longitudinal force exerted on the wheel at the contact point N 
F, vertical load on wheel N 


2.2. Electric motorcycle (EM) 


The required components to develop EM are quite similar to ICEM by which the electric motorcycle 
or EM also has a nonlinear rotational damper, inertia block, differential, vehicle body, and tires. However, 
there are slight differences in developing the EM model where the additional components are the controlled 
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PWM voltage, H-bridge, DC motor, and gearbox. In addition, the generic engine is replaced by a generic 
battery. The power to run the EM is produced by the generic battery where batteries are the fundamental 
components in every electronic application. The equation for the generic battery is shown as (24). The 
nominal voltage and ampere-hours (AH) rating input can be varied. 


yaw -G2 o» 


Next, the controlled PWM voltage will receive the voltage from the generic battery. The PWM 
voltage is determined based on the reference voltage. As a result, the demanded duty cycle can be calculated. 
The equation for the demanded duty cycle is shown as (25). Then, controlled PWM voltage is connected with 
the H-bridge. The H-bridge is used as a motor driver. The H-bridge will accept the PWM voltage as an input. 
The equation for the H-bridge is shown as (26). 


S Vref -Vmin 
D = 100 x (Fermin) % 2) 
V-Vpwm 
H, = A P - loutRon (26) 
pwm 


The DC motor is connected to the output of the H-Bridge. The DC motor is modeled with two 
characteristics namely electrical and torque. The equations for DC motor are shown as (27)-(31). 


E 27) 
Te = k¿i (28) 
Tew = vi (29) 
k,iw = kwi (30) 
ky =k, (31) 


Lastly, the gearbox functions as a gear ratio. The gear ratio is the ratio between the input shaft and 
output shaft angular velocity. The input and output connections of the gearbox are connected to the inertia 
block in order to ensure the motorcycle accelerates in a positive direction. The equations for the gearbox are 
shown as (32)-(35). Table 1 summarizes the parameters used for modeling the EM. 


w, =N-w, (32) 
T2 =N-T, (33) 
P =w: Ti (34) 
P, = —w, Ty (35) 


3. SIMULATION DESIGN 

All the equations in section 2 are converted into a Simulink block diagram. The type of generic engine 
used in the ICEM system is a spark-ignition engine whereby its power demand coefficient (p1, p2, and p3), is 
set by 1. In the subsystem, the gear ratio can be adjusted by varying the number of teeth of the base on gear, NB 
and number of teeth of follower on gear, NF. Furthermore, the rotation of the output shaft can be selected either 
in the opposite or in a similar direction towards the input shaft by using the simple gear. For this system, a 
selection is made whereby rotation of the output shaft moves in the same direction towards the input shaft. 

In designing the EM system, it is similar to the ICEM system with the differences in the system 
design of EM as compared to ICEM is the generic battery, controlled PWM voltage, H-bridge, DC motor and 
gearbox subsystem. In the EM system is a generic battery which is significant in supplying battery nominal 
voltage, Vo, and is controlled by PWM voltage to loads likes DC motor based on the reference voltage. The 
EM system includes the gearbox to allow a variety of gear ratio combinations. The input and output 
parameters of the models are shown in Table 2. 
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Table 2. Parameters of simulation for ICEM and EM system 


Symbol Parameter Value and unit 
g gravitational acceleration 9.81 m/s” 
B inclination angle 0° 
Vw headwind speed 0 km/hr 
m ICEM mass 165kg 
EM mass 80kg 
h height of vehicle CG above the ground 0.5m 
a distance of front axles, respectively, from the normal 1.4m 
projection point of vehicle CG in the common axle 
b distance of rear axles, respectively, from the normal 1.6m 
projection point of vehicle CG in the common axle 
n number of wheels on each axle 1 
A effective frontal vehicle cross-sectional area 3m? 
Ci aerodynamic drag coefficient 0.4 
p mass density of air 1.18kg/m° 
V, nominal voltage 72V 
x ratio of ampere-hours (AH) left to the number of AH 60AH 
Vinin minimum reference voltage OV 
Vine maximum reference voltage 72V 
k, back emf constant 7.2e-5150 
N gear ratio 10 
bi, b2, ...b5 damping coefficient le-4, 0, 1e-5, 0, 1e-6 
Pore maximum power from the engine 1.0e+5W 
Q engine speed 800rpm 
P1P2P3 power demand coefficient for spark-ignition engine 1,1,1 
SEB follower-base gear ratio 2 
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4. RESULT AND DISCUSSION 

In the simulation, both motorcycle models are exposed to different road conditions namely, dry, and 
wet, to investigate the effect of these conditions on to slip ratio of the motorcycle. Figure 2 (a) illustrates the 
results of the ICEM and EM front and rear wheel slip ratio on a dry road. From this figure, for the front 
wheel of ICEM, the slip ratio is initially recorded at 62.83 before slowly declining and remains constant at a 
slip ratio of 0.035. For the rear wheel slip ratio of ICEM, the slip ratio is recorded to be 16.87 initially before 
finally remaining constant with a slip ratio of 0.002. 

While for the EM on the dry road, on the front wheel the slip ratio at the beginning is recorded at 0 
before a slight increase to 0.27 at 0.02 seconds and gradually decreases and remains constant with 0.058. For the 
rear wheel slip ratio of EM, the slip ratio only shows a slight increase from 0 initially to finally remain constant 
at 0.00167. Both the ICEM and EM are then simulated under the wet road condition and are documented in 
Figure 2 (b). Based on Figure 2 (b), the slip ratio of the front wheel of ICEM is sharply declined when initially 
peaked at 62.83 before finally remaining constant at 0.157. Whereas the final value of the ICEM’s rear-wheel 
slip ratio is at 0.072 when initially is recorded at 21.07 during the beginning of the simulation. 
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Figure 2. ICEM and EM front and rear wheel slip ratio on (a) dry road and (b) wet road 


On the other hand, the front wheel slip ratio of EM is recorded at 0 initially with a slight increase in 
the value of 3.93 between 0 and 0.07 seconds before gradually decreases and remain stable at a ratio of 
0.304. For the rear wheel slip ratio of EM, the slip ratio started at a ratio of 0 and a slight increase to 0.64 at 
1.23 seconds before gradually declined and settled at a slip ratio of 0.2. From the Figures 2 (a) and 2 (b), the 
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slip ratio for ICEM on the front and rear wheel for both road conditions shows that the slip ratio for ICEM is 
higher than EM. 

The speed profiles of both ICEM and EM under a dry road condition are shown in Figure 3 (a) 
whereas Figure 3 (b) illustrates the speed profile of both systems when subjected to wet and dry road 
conditions. In terms of the vehicle's speed, under the dry road condition, the ICEM is able to reach up to 
91.32 km/h while the EM managed to rise to the speed of 74.77 km/h within the 20 seconds of the simulation 
time. When the systems are subjected to the wet road condition, the speed of the ICEM is observed to reach 
to 82.55 km/h while the speed of the EM is noted to reach 64.96 km/h. 

Note that the speed of the ICEM is higher than the speed of the EM since the slip ratio is inversely 
proportional to the friction coefficient and directly proportional to motorcycle speed. The higher the slip 
ratio, the higher the motorcycle speed. Furthermore, the graph pattern for slip ratio of both ICEM and EM on 
dry and wet road conditions indicates a similar pattern as through time, the slip ratio slowly remains stable. 
As the slip ratio is stable, the motorcycle will reach a steady condition and a suitable motorcycle speed when 
driven on dry and wet road conditions. 


Speed (km/h) 


o 2 4 6 8 10 12 14 16 18 20 o 2 4 6 8 10 12 14 16 18 20 
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(a) (b) 
Figure 3. ICEM and EM speed on (a) dry road and (b) wet road 
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Figure 4. Relationship between slip ratio and speed of ICEM through the time on (a) dry road and 
(b) wet road 


To acquire an understanding in regard to the relationship between slip ratio and motorcycle's speed 
through the times, the three-dimensional (3D) graphs are being shown in Figures 4 (a), 4 (b), and 5 (a), 5 (b). 
Through all those 3D figures, it can be observed that within the same time frame the slip ratio of both ICEM and 
EM reach a constant value. The slip ratios are recorded to increase stably until the maximum speed is reached. 
Based on these findings, in the dry road condition, the maximum speed the ICEM can travel is around 91.32 
km/h while the maximum speed of EM is approximately 74.77 km/h. In the wet condition, the maximum speed 
is less for both motorcycles with the maximum speed for ICEM is 82.55 km/h and EM is 64.96 km/h. 

From Figure 4 (a), at the maximum speed on the dry roadway, the slip ratio is recorded at 0.035 
while in Figure 4 (b) at the maximum speed on the wet roadway, the slip ratio is observed at 0.157. A higher 
slip ratio on the wet roadway is expected as the tire friction reduces in this road condition. The same evidence 
can be seen in the case of EM. In Figure 5 (a), at the maximum speed on the dry road the slip ratio is 
recorded at 0.058 whereas in Figure 5 (b) at the maximum speed on the wet roadway, the slip ratio is 
recorded at 0.304, a higher value than the slip ratio on a dry road. It is expected that the slip ratio of the 
vehicle is much higher on the wet roadway as the tire friction reduces in this road condition. 
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Figure 5. Relationship between slip ratio and speed of EM through the time on (a) dry road and (b) wet road 


5. CONCLUSION 

There are numerous studies have been attentively reviewed and studied to provide an explicit 
explanation in regards to the suitable speed range of ICEM and EM when driven on dry and wet road 
conditions. From this study, it can be concluded that the maximum speed range for ICEM to travel on a dry 
road condition is up to 91.32 km/h whereby the calculated speed for wet road conditions is 82.55 km/h. 
Meanwhile, for EM, the maximum speed when driven on a dry road condition is found to be 74.77 km/h as 
compared to when driven on a wet road condition which is 64.96 km/h. 
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